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ABSTRACT 

The mass function and statistics of binaries provide important diagnostics of the star formation process. 
Despite this importance, the mass function at low masses remains poorly known due to observational difficulties 
caused by the faintness of the objects. Here we report the microlensing discovery and characterization of 
a binary lens composed of very low-mass stars just above the hydrogen-burning limit. From the combined 
measurements of the Einstein radius and microlens parallax, we measure the masses of the binary components 
of 0. 10 ± 0.01 Mq and 0.09 ± 0.01 M Q . This discovery demonstrates that microlensing will provide a method 
to measure the mass function of all Galactic populations of very low mass binaries that is independent of the 
biases caused by the luminosity of the population. 



1. INTRODUCTION 
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2. INTRODUCTION 

Microlensing occurs when a foreground astronomical ob- 
ject (lens) is closely aligned to a background star (source) and 
the light from the sour ce star is deflected by the gravity of 
the lens (Einstelnl ll936l) . The phenomenon causes splits and 
distortions of the source star image. For source stars located 
in the Milky Way Galaxy, the separation between the split im- 
ages is of order milli-arcsecond and thus the individual images 
cannot be directly observed. However, the phenomenon can 
be photometrically observed through the brightness change 
of the source star caused by the ch ange of the relative lens- 
source s eparatio n (Paczvriski 1986). Since the first discovery 
in 1993 (lAlcock et alJl!993b Iudalskilll993l) . there have been 
numerous detections of microlensing events toward the Large 
and Small Magellanic Clouds, M31, and mostly the Galac- 
tic bulge fields. Currently, microlensing event s are being de- 
tected at a rate o f nearly 1000 events per year (Udalskil l2008l ; 
Bon d et alJl200l . 

The properties of multiple systems such as binary frequency 
and mass function provide important constraints for star for- 
mation theories, enabling a concrete, qualitative picture of the 
birth and evolution of stars. At very low masses down to and 
below the hydrogen burning minimum mass, however, our 
understanding of formation processes is not clear due to the 
difficulties of observing these objects. Over the last decade, 
there have been se veral s earch e s for very low mass bina - 
ries [see reviews b yjBasri (2000); Oppenheimer et al. (2000); 
lKirkpatrici3d2005b ; lB"urgasser et alj d2Q07l) 1. Despite these ef- 
forts, the number of very low mass binarie^l is not big enough 
to strongly constrain their formation processes. 

Microlensing occurs regardless of the brightness of lens ob- 
jects and thus it is potentially an effective method to inves- 
tigate the mass function of low-mass binaries. For lensing 
events caused by single-mass objects, it is difficult to mea- 
sure the lens mass because the Einstein time scale t^, which is 
the only observable that provides information about the lens 
for general lensing events, results from the combination of 
the mass and distance to the lens and the transverse speed 
between the lens and source. The degeneracy can be par- 
tially lifted by measuring either an Einstein radius or a lens 
parallax and can be completely broken by measuring both. 
Einstein radii are measured from the deviati on in lensin g 
light curves caused by the finite-source effect (iGouldll 19941) . 
Most microlensing events produced by binaries are identi- 
fied from the anomalies involved with caustic approaches 
or cr ossings during which the finite-source effect is impor - 
tant dNemiroff & Wickramasinghe1ll994l; IWitt & Maolll994l) . 
Therefore, Einstein radii can be routinely measured for the 
majority of binary-lens events. The microlens parallax is de- 
fined by 
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where 7r re i = AU(Dj" 1 -D$ l ) is the lens-source relative parallax 
and Dl and Ds are the distances to the lens and source, re- 
spectively. In general, parallaxes can be measured for events 
that last long enough that the Earth's motion can no longe r 
be approximated as rectilinear during the event (IGouldll 9921) . 
The chance to measure the lens parallax for binary-lens events 
is higher than that of single-lens events because the average 
mass of binaries is bigger and thus time scales tend to be 
longer. In addition, the well-resolved caustic-crossing part of 
lensing light curves provide strong constraints on the lensing 
parameters and thus helps to pin-down enough anchor points 
on th e light curve to extr act otherwise too-subtle parallax ef- 
fects (lAn & GoulJ200l) . 

The number of binary-lens events with well-resolved 
anomalies is increasing with the advance of observational 
strategies such as the alert system and follow-up observa- 
tions. The increase of the monitoring cadence of existing and 
planned survey experiments will make the detection rate even 
higher. Although binary microlensing is biased toward sepa- 
rations similar to the Einstein radius, it is easily quantifiable 
for next-generation experiments that have continuous "blind" 
monitoring. Therefore, microlensing will be able to provide 
an important method to discover very low mass binaries and 
to investigate their mass function. 

In this paper, we present the microlensing discovery and 
characterization of a very low mass binary. We use this 
discovery to demonstrate that microlensing will provide a 
method to measure the mass function of very low-mass bina- 
ries that is free from the biases and difficulties of traditional 
methods. 

3. OBSERVATIONS 

The low-mass binary was detected from observations of 
the microlensing event OGLE-2005-BLG-153. The event 
occurred on a Galactic bulge star located at right ascen- 
sion a = 18 h 04 m 17 s .30 and declination S = -28°40'49".3 
(J2000), which corresponds to the Galactic coordinates (l,b) = 
(2.3°, -3.4°). The event was detected by the Optical Gravita- 
tional Lensing Experiment (OGLE) using the 1.3 m Warsaw 
telescope of Las Campanas Observatory in Chile and was an- 
nounced as a probable microlensing event on 14 April 2005. 
The event was independently identified by the Microlensing 
Observation in Astrophysics (MOA) as MOA-2005-BLG-023 
using the 0.6 m of Mt. John Observatory in New Zealand. 

An anomaly alert was issued on 3 July 2005 by the Probing 
Lensing Anomalies Network (PLANET) collaboration. Fol- 
lowing the alert, the PLANET, RoboNet, and Microlensing 
Follow-Up Network (//FUN) teams intensively observed the 
event by using eight telescopes located on three different con- 
tinents. These telescopes include the PLANET 1.54 m Danish 
Telescope of La Silla Observatory in Chile, PLANET 1.0 m 
of Mt. Canopus Observatory in Australia, PLANET 0.6 m of 
Perth Observatory in Australia, //FUN 1.3 m SMARTS tele- 
scope of CTIO in Chile, RoboNet 2.0 m Faulkes Telescope 
S. (FTS) in Australia, RoboNet 2.0 m Faulkes Telescope N. 
(FTN) in Hawaii, and RoboNet 2.0 m Liverpool Telescope 
(LT) in La Palma, Spain. Thanks to the follow-up observa- 
tions, the light curve was densely resol ved. 

The event was analyzed before by iSkowron et al.l ( 120071) 
based on the data from OGLE observations as one of 9 binary- 
lens events detected in 2005 season. Here we reanalyze the 
event in depth with the addition of data from follow-up obser- 
vations focusing on the physical parameters of the lens sys- 
tem. 
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FIG. 1 . — Light curve of the microlensing event OGLE-2005-BLG-153. Also presented is the model curve for the best-fit solution. The upper panel shows a 
zoom of the caustic-induced perturbation region. 



4. MODELING 

Figure 1 shows the light curve of the event. It is character- 
ized by the sharp rise and fall occurring at around the helio- 
centric Julian dates (HJD) of 2453556 and 2453560. These 
features are caused by the crossings of the source star across 
a caustic, which represents a set of source positions at which 
the lensing magnification of a point source becomes infinite. 
Therefore, the existence of such a feature immedia tely reveals 
that t he lens is composed of two masses dMao & Paczynskil 
fl99l . 

Characterization of binary lenses requires modeling of lens- 
ing light curves. We test 3 different models. In the first model, 
we test a static binary model (standard model). In this model, 
the light curve is characterized by 7 parameters. The first set 
of three parameters are needed to describe the light curves of 
single-lens events: the time required for the source to transit 
the Einstein radius, ?e (Einstein time scale), the time of the 
closest lens-source approach, to, and the lens-source separa- 
tion in units of the Einstein radius at the time of to, uo (impact 
parameter). Another set of three parameters are needed to de- 
scribe the deviation caused by the lens binarity: the projected 
binary separation in units of the Einstein radius, s, the mass 
ratio between the binary components, q, and the angle of the 
source trajectory with respect to the binary axis, a (source tra- 
jectory angle). Finally, an additional parameter of the ratio of 
the source radius to the Einstein radius, p+ = 6+/6-E (normal- 
ized source radius), is needed to incorporate the deviation of 
the light curve caused by the finite-source effect. In the sec- 



ond model, we consider the parallax effect by including two 
additional parallax parameters of ttb^i and 7r E £, which are the 
two components of the microlensing parallax vector 7Te pro- 
jected on the sky in the direction of north and east celestial 
coordinates. In the last model, we additionally check the pos- 
sibility of the effect on the lensing light curve caused by the 
orbital motion of the lens. The orbital motion affects the lens- 
ing magnifications in two different ways. First, it causes the 
binary axis to rotate or, equivalently, makes the source trajec- 
tory angle change in time. Second, it causes the separation 
between the lens components to change in time. The latter ef- 
fect causes alteration of the caustic shape in the course of an 
event. To first order, the orbital effect is parameterized by 

a(t) = a(to)+u(t^) (2) 

and 

s(t) = s(to)+s(*-^) (3) 

where the orbital parameters lu and s represent the rates of 
change of the source trajectory angle and the projected binary 
separation, respectively. Considering the orbital effect is im- 
portant not simply to constrain the orbital motion of the lens 
system but also to precisely determine the lens mass. This 
is because both the motions of the observer (parallax effect) 
and the lens (orbital effect) have a similar effect of causing 
deviations of the source trajectory from a straight line. Then, 
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TABLE 1 
Fit Parameters 





standard 


parallax 




parallax+orbit 






(«o > 0) 


(H < 0) 


(un > 0) 


(un < 0) 


2 

X 


^ in/ 010 

2406.812 


1836.199 


1982.545 


1672.572 


1674.966 


to (HJD') 


3549.8286±0.0556 


3551.2263±0.0950 


3549.4179±0.0692 


3549.1650±0.0992 


3548.6344±0.1452 


"0 


0.502±0.002 


0.481 ±0.003 


-0.515±0.002 


0.566±0.003 


-0.572±0.004 


t E (days) 


46.428±0.lll 


48.735±0.149 


50.019±0.229 


44.541±0.121 


44.952±0.173 


s 


0.837±0.00l 


0.836±0.001 


0.832±0.001 


0.848±0.001 


0.847±0.001 




0.848±0.003 


0.848±0.003 


0.833±0.003 


0.871±0.003 


0.865±0.002 


a 


5.902±0.00l 


5.956±0.003 


O.373±0.003 


5.900±0.002 


0.396±0.003 


P* 


0.0l7±0.00l 


0.017±0.001 


0.016±0.001 


0.018±0.001 


0.018±0.001 


Vi 


0.57l±O.Ol9 


0.510±0.016 


0.545±0.016 


0.501±0.014 


0.495±0.013 






0.762±0.040 


O.555±0.078 


0.106±0.047 


0.196±0.175 






0.699±0.031 


0.921±0.042 


0.419±0.031 


0.385±0.029 


S 








-0.0651±0.0020 


-0.0693±0.0022 


a 








0.0064±0.0068 


-0.0004±0.()101 



NOTE. — HJD' = HJD -2450000. The parameters of the best-fit model are marked in bold font. 



if the orbital motion of the lens is not considered despite its 
non-negligible effect, the deviation of the lensing light curve 
caused by the orbital effect may be explained by the paral- 
lax effect. This will cause a wrong determination of the lens 
parallax and the resulting lens mass. 

When either the effect of the parallax or orbital motion is 
considered, a pair of source trajectories with impact parame- 
ters uo > and uq < result in slightly different light curves 
due to the breakdown of the mirror-image symmetry of the 
source trajectory with respect to the binary axis. We, there- 
fore, check both models with uo > and uo < whenever the 
parallax or orbital effect is considered. As a result, the total 
number of tested models is 5. 

To find the best-fit solution of the lensing parameters, we 
use a combination of grid and downhill approaches. It is 
difficult to find solutions from pure brute-force searches be- 
cause of the sheer size of the parameter space. It is also diffi- 
cult to search for solutions from a simple downhill approach 
because the \ 2 surface is very complex and thus even if a 
solution that apparently describes an observed light curve is 
found, it is hard to be sure that all possible \ 2 minima have 
been searched. To avoid these difficulties, we use a hybrid 
approach in which a grid search is conducted over the space 
of a subset of parameters (grid parameters) and the remain- 
ing parameters are searched by a down-hill approach to yield 
minimum \ 2 at each grid point. Then, the best-fit solution 
is found by comparing the \ 2 values of the individual grid 
points. We set s, q, and a as grid parameters because they are 
related to the features of light curves in a complicated pattern 
while other parameters are more directly related to the identi- 
fiable light curve features. For the down-hill \ 2 minimization, 
we use a Markov Chain Monte Carlo method. 

5. RESULTS 

In Table 1, we present the results of modeling along with 
the best-fit parameters for the individual models. It is found 
that the effects of parallax and orbital motion are needed to 
precisely describe the light curve. We find that the model 
with the parallax effect improves the fit by A\ 2 = 571. The 
fit further improves by A\ 2 =164 with the addition of the or- 
bital effect. We note that the values of the parallax parameters 
from the "parallax+orbit" model are different from those de- 
termined from the "parallax" model. This demonstrates that 
consideration of the orbital effect is important for the precise 
measurement of the lens parallax. 
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FIG. 2. — Geometry of the binary lens system responsible for the lensing 
event OGLE-2005-BLG-153. In the lower panel, the two filled dots repre- 
sent the locations of the binary lens components. The dashed circle is the 
Einstein ring corresponding to the total mass of the binary. The ring is cen- 
tered at the center of mass of the binary (marked by '+'). The line with an 
arrow represents the source trajectory. We note that the trajectory is curved 
due to the combination of the effects of parallax and lens orbital motion. 
The closed figure composed of concave curves represents the positions of the 
caustic formed by the binary lens. All lengths are normalized by the Einstein 
radius. The temperature scale represents the magnification where brighter 
tones imply higher magnifications. The upper panel shows a zoom of the 
boxed region. We note that the caustic shape slightly changes due to the or- 
bital motion of the lens. We present the caustics at the two different moments 
of the caustic entrance and exit of the source star. 

In Figure [1] we present the model light curve on the top 
of observed data points. Figure [2] shows the geometry of the 
lens system corresponding to the best-fit solution, i.e. "paral- 
lax+orbit" model with uq > 0. In the figure, the filled dots 
represent the locations of the binary components, the dashed 
circle is the Einstein ring corresponding to the total mass of 
the lens, the closed figure composed of concave curves is the 
caustic formed by the lens, and the curve with an arrow rep- 
resents the source trajectory. The upper panel shows an en- 
largement of the region where the source trajectory crosses 
the caustic. The shape of the caustic changes in time due to 
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FIG. 3. — Positions of the source (lensed) star with respect to the centroid 
of clump giants in the instrumental color-magnitude diagram. 

the orbital motion of the lens, and thus we present the caustics 
at two different moments of the caustic entrance and exit of 
the source star. 

Among the two quantities needed for the determination of 
the lens mass, the microlens parallax is obtained directly from 
the parallax parameters determined from modeling by 



7T E = (-7T] 



2 

EJV 



+ 7T 



k V2 



E.EJ 



= 0.432 ±0.042. 



(4) 



On the other hand, the Einstein radius is not directly obtained 
from modeling. Instead, it is inferred from the normalized 
source radius p+, which is determined from modeling, com- 
bined with the information about the angular source radius 8+. 
The angular source radius is determined from the information 
of the de-reddened color of the source star measured by us- 
ing the centroid of clump giant stars in the color-magnitude 
diagram as a reference position under the assumption that the 
source a nd clump centroi d experience the same amount of ex- 
tinction dYoo et al.l 120041) . Figure [3] shows the instrumental 
color-magnitude diagram constructed based on CTIO V and / 
band images and the locations of the source star and the cen- 
troid of clump giants. By measuring the offsets in the color 
and magnitude between the source and the centroid of clump 
giants combined with the known color and absolute magni- 
tude of the clump centroid of [(V-7),Af/] c = (1.04,-0.25), 
we estimate that the de-reddened magnitude and color of the 
source star are 1q = 13.16 and (V— I)o= 1 .09, respectively, im- 
plying that the source is a clump giant with an angular radius 
of 9+ = 11.72 ± 1.01 /xas. Here we adopt a Galactocentric 
distance of 8 kpc and the offset of the bar of the field is 0.4 
kpc toward the Sun and thus the distance to t he clump cen- 
troid i s 7.6 kpc based on the Galactic model of H an & Gould! 
( 12003b . Then, with the measured normalized source radius of 
Pi, = 0.018 ±0.001, the Einstein radius is estimated as 



— =0.66 ±0.06 mas. 



(5) 



Together with the Einstein time scale, the relative lens-source 
proper motion is obtained by p = #e Ae = 5 .38 ± 0.47 mas yr" 1 . 
With the measured Einstein radius and lens parallax, the 



K1T E 



:0.19±0.02M, 



0> 



(6) 



where k = 4G/(c 2 AU). With the known mass ratio between 
the binary components, the masses of the individual binary 
components are determined, respectively, as 



1 



and 



m 2 = 



l+q 



-m = 0.10 ± 0.01 Mr, 



© 



l+q 



m = 0.09 ±0.01 M, 



©• 



(7) 



(8) 



This implies that both lens components are very low mass 
stars with masses just above the hydrogen-burning limit of 
0.08 Mq. The distance to the lens is determined as 

D L =— ^ — = 2.42 ±0.21 kpc, (9) 

where its = AU / Ds is the parallax of the source star. From this 
distance to the lens combined with the Einstein radius, it is 
found that the two low-mass binary components are separated 
with a projected separation of 



r± = sDl9e = 1.35 ±0.12 AU. 



(10) 



It is also found that the lens velocity in the frame of the lo- 
cal standard of rest is v = (v±,V||) = (-20.9 ± 31.9, 15.1 ± 
31.9) km s" 1 , where vj_ and vm are the velocity components 
normal to and along the Galactic plane, respectively. We note 
that the errors in v are dominated by the unknown proper mo- 
tion of the source, which is assumed to be ± 100 km s" 1 in 
the Galactic frame. The velocity and the distance to the lens 
imply that the lens is in the Galactic disk. 

In addition to the parallax effect, the relative lens-source 
motion can, in principle, also be aff ected by the orbital motion 
of the source star if it is a binary ( Smi th. Mao. & Paczynskil 
2003). We check the possibility that this so-called "xallarap" 
(reverse of "parallax") effect influences the parallax determi- 
nation. For this, we conduct additional modeling including 
the xallarap effect. For the description of the xallarap effect, 
3 additional parameters of the phase angle, inclination, and or- 
bital period are needed under the assumption that the source 
moves in a circular orbit. From this analysis, we find that 
the xallarap effect does not provide a better model than the 
parallax model. In addition, the best-fit occurs for an orbital 
period of ~ 1 yr, which corresponds to the orbital period of 
the Earth around the Sun. Furthermore, the best-fit values of 
the inclination and the phase angle are similar to the eclip- 
tic longitude and latitude of the source star. All these facts 
imply that the xallarap interpretation of the light-curve de- 
viation is less likely and support the pa rallax interpretation 
dPoindexter et alJl2005tlDonp et al.ll2009h . 

From the orbital parameters uj and s determined from mod- 
eling along with the assumption of a circular orbit, one can 
obtain the usual orbital parameters of the semi-major axis, a, 
orbital period, P, and inclination, i, of the orbit of the binary 
lens from the relations 



P = 



1/2 



cos i = -(xB) 



1/2 



(ID 



Here B = r\u 2 /{Gm), x = sin (f), and is the angle between 
the vector connecting the binary components and the line of 
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sight to the lens such that the projected binary separation is 
r± = a sin <fi. The value of x is obtained by solving the equation 
x 3 -B x 2 -x + (A 2 + 1)B = 0, where A = {s/s)/uj ( iDong et all 
120091) . We find that the semi-major axis is a = 1 .46 ± 0.08 AU 
and the period is P = 4.05 ± 0.19 yrs. The inclination of the 
orbital plane is i = 88.3° ± 1.1°, implying that the orbit is very 
close to edge on. 

We can also constrain the surface brightness profile of the 
source star by analyzing the caustic-crossing parts of the light 
curve. We model the source brightness profile as 



F 



i-r> 



, 3 

1 COS I 

2 



(12) 



where T\ is the linear limb-darkening coefficient and 9 is the 
angle between the normal to the stellar surface and the line 
of sight toward the source star, and F is the source flux. We 
measure the coefficient in / band of T/ = 0.501 ± 0.014. The 
measured coeffi cient is cons istent with the theoretical value 
of clump giants (IClaretll2000t) . 

6. DISCUSSION AND CONCLUSION 

We analyze the light curve of a binary-lens microlensing 
event OGLE-2005-BLG-153, which exhibits a strong caustic- 
crossing structure on the light curve. By measuring both the 
Einstein radius and the lens parallax, we could uniquely mea- 
sure the masses of the individual lens components. The mea- 
sured masses were 0.10 ±0.01 M and 0.09 ±0.01 M , re- 
spectively, and thus the binary was composed of very low- 
mass stars just above the hydrogen-burning limit. 

Although the event OGLE-2005-BLG-153 is one of few 
cases with well-measured lens masses among the 5000 mi- 
crolensing events discovered to date, the event characteris- 
tics that enabled this mass measurement are likely to become 
common as next-generation microlensing experiments come 



on line. Because next-generation experiments will provide in- 
tense coverage from sites on several continents, most caustic- 
crossing binaries will yield masses. Moreover, because next- 
generation cadences will be independent of human interven- 
tion, rigorous characterization of the selection function will 
be straightforward. Finally, for reasonable extrapolations of 
the mass function of stars close to and below the hydrogen- 
burning limit, we can anticipate an important fraction of the 
roughly thousand events per year expected to be detected from 
next-generation surv eys to be d ue to low-mass objects includ- 
ing brown dwarfs (Gould 2009). Hence, the mass function, 
at least of objects within binaries, will be measurable for all 
Galactic populations of low-mass stellar and substellar objects 
in the near future, independent of biases caused by the lumi- 
nosity of the population. 
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